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^ . Abstract 

We investigate some connections between the continuum and atomistic descriptions of de- 
formable crystals, using some interesting results from number theory. The energy of a deformed 
crystal is calculated in the context of a lattice model with binary interactions in two dimensions. 
A new bond counting approach is used, which reduces the problem to the lattice point problem 
of number theory. When the crystal shape is a lattice polygon, we show that the energy equals 
the bulk elastic energy, plus the boundary integral of a surface energy density, plus the sum over 
the vertices of a corner energy function. This is an exact result when the interatomic potential has 
finite range; for infinite-range potentials it is asymptotically valid as the lattice parameter zero. 
The surface energy density is obtained explicitly as a function of the deformation gradient and 
boundary normal. The corner energy is found as an explicit function of the deformation gradient 



and the normals of the two facets meeting at the corner. For more general convex domains with 
possibly curved boundary, the surface energy density depends on the unit normal in a striking 
way. It is continuous at irrational directions, discontinuous at rational ones and nowhere differ- 
entiable. This pathology is alarming since it renders the surface energy minimization problem 
(under domain variations) ill-posed. An alternative approach of defining the continuum region 
$h \ is introduced, that restores continuity of the surface energy density function. 



1 Introduction 

This article is concerned with the derivation of continuum surface energy from a standard lattice 
model, by exploiting results related to certain lattice point problems of number theory, e.g. ||BU 
IMllHullIKMl giTl. 

The problem of relating discrete energies of crystals to their continuum counterparts has a much 
simpler analogue: that of mass. Suppose the crystal is two-dimensional and consists of those atoms 
(points in the lattice L = Z 2 ) that belong to a compact set f2. Consider the problem of measuring the 
mass m of the crystal, given that each atom has unit mass. Then the exact mass is m = fl L). 
Suppose further that the crystal is faceted; in particular, let Vt be a lattice polygon: one whose vertices 
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are lattice points. This implies that the facets, or sides of are along crystallographic directions. 
One would be tempted to write the mass asm = J n ldx = |0| independently of f2, but it turns out 
that this is not true in general, though it holds in the limit as the crystal becomes large. The exact 
mass can be found using a remarkable and quite old result, namely Pick's Theorem JH1 (1899): 
The area of a lattice polygon equals the number of atoms in the interior, plus half the number of 
atoms on the boundary, minus one. This provides a discrete expression for the area — a "continuum" 
quantity — but it is equivalent to the following continuum expression (see Lemma l2Tl below) for the 
number of atoms (a "discrete" quantity) or the mass: 



where n is an outward normal to dVt whose components on each facet are irreducible integers, 
and Oi are the dihedral angles of the N vertices of Vt. This is a continuum expression of the mass 
as a bulk integral of a (unit) mass density, plus the boundary integral of a surface mass density, 
plus contributions from the corners. In this paper we study the energy of a crystal due to pairwise 
atomic interactions, with emphasis on its dependence on the shape of the boundary. Some of our 
results express this energy of faceted crystals in a similar, continuum form, as bulk energy, plus the 
boundary integral of a surface energy density, plus a sum of corner energies. The surface and corner 
energies are found explicitly in terms of the interatomic potential, the deformation gradient, and the 
normal vectors. We also study the surface energy for smooth boundaries. 

Perhaps the first rigorous derivation of continuum energy functions from atomistic models is 
due to Blanc, Le Bris and Lions [BBLJ, who study (among other problems) the energy of a crystal 
governed by an interatomic potential with binary central-force interactions, subject to a prescribed 
smooth deformation. An asymptotic form of the energy is obtained as the lattice parameter e tends 
to zero. The dominant term is the usual elastic energy J n W(Wy(x) )dx, where VL is the macroscopic 
reference domain occupied by the crystal, y : — > M 3 is the prescribed smooth deformation, and 
the elastic stored energy function W is determined by the interatomic potential (cf. (13.31) below). 
The next term, of order e in Theorem 3 of [BBLJ, is a surface integral over OVl that involves values 
of the deformation gradient and unit normal. The form of this surface energy is not explicit and 
it is not clear to what extent it can be expressed as a function of those variables. Terms of order 
e 2 include a volume integral of higher gradient energy, but also surface terms; the latter are left 
unspecified. 

In studying the energy, one encounters sums of the form e 3 XLe(nneL) f( x )- Here eL is the 
rescaled lattice L = Z 3 with lattice parameter e. While these Riemann sums converge to f Q f as 
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e — > 0, there are higher order terms. Often, these are surface terms of order e, the same order as 
the surface energy in HBBLL The asymptotic behavior of Riemann sums, even in simple smooth 
domains, is not completely understood at present, unless strong specializing assumptions are made 
[BRJ. This is true even when / is constant, and the problem reduces to the lattice point problem, 
well-known in number theory: find the asymptotics of the number of lattice points in Q, namely 
fl eL) as e — > 0, equivalently #(rfi fl L) as r — > oo [BRl lIKMllTsll . Even in two dimensions, 
in the seemingly simple case where f2 is the unit disk, the problem — initially posed by Gauss — 
is not completely settled. Currently, the best available estimate seems to be due to Huxley HHuL 
Letting R(r) = #(rf2 PI L) — \rVt\ be the remainder, then for the disk, or for sets with positive 
boundary curvature, R(r) = 0(r a (\ogr) b ), with a = §§ « 0.63 and b = « 2.25 tiHul The 
remainder depends strongly on the shape of <9f2; for f2 a lattice polygon (one whose vertices are 
lattice points), R(r) = 0(r). Letting e = 1/r, the remainder would corresponds to a higher order 
term Q(e) = e 2 R(l/e) in the Riemann sum (of a constant function) e 2 ^(VtneL) = +Q(e). For 
smooth, strictly convex domains, the order of Q(s) is between e and e 2 [TKMJ, that is, between the 
orders of the surface and the gradient energy of [BBLJ. In addition, Q(e) is highly oscillatory and 
very difficult to characterize explicitly. In [BBLJ it is assumed that there is a sequence e = — > 
as the integer k — > oo, such that the number of lattice points in f2 is precisely (in three 

dimensions). This means that for some sequence — > oo, #(rfefi D L) = |rfeQ|. This is equivalent 
to Q(Ek) = 0, thus it eliminates certain undesirable higher order terms from a Riemann sum of 
the elastic energy. Unfortunately however, it is not known for which choices of domain such a 
sequence exists. As shown in one dimension by Mora-Corral HMoll , the higher order terms in the 
asymptotic expansion of the energy depend on the choice of the sequence e^. 

For lattice polygons, Q(e) = 0(e) and it can be characterized explicitly in terms of geometrical 
quantities of f2. The dominant term of Q(e) is a surface integral over dVt (the second term in (11.11) 
after scaling); as such it furnishes an explicit contribution to surface energy that plays a rather subtle 
role, as we show in the present paper. 

Crystals often occur in faceted form in their natural state. This means that they are polyhedral, 
and the facets inhabit crystallographic planes. This is because of surface energetics affecting crystal 
growth, but also because of cleavage fracture that creates new surfaces along special lattice planes. 
It turns out that considerable progress has been made in the lattice point problem for regions that are 
lattice polytopes, i.e., polyhedra whose vertices are lattice points HBPllBRllRell . For lattice polygons 
in two dimensions the problem was solved by Pick [PiJ, however, in three or more dimensions 
some issues remain open [BPl IBRI . In Section [2] we assume that Vt is a lattice polytope. In order 
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to isolate the surface energy from other effects, such as higher gradient energy, we assume that the 
deformation is homogeneous (y is affine). To keep the geometry simple, we confine our analysis to 
two dimensions. Unlike [BBLl iMol . initially we do not employ a limit process, but rather a bond 
counting technique. For each lattice vector w E L, we find the number of bonds between points in 
OnL that are w apart. This calculation reduces to a number of lattice point problems. 

In Section [3] we compute the energy of polygonal (faceted) crystals. For an interatomic potential 
of finite, but otherwise arbitrary range, we compute the energy of essentially any convex lattice poly- 
gon exactly (Proposition 13.11) . This result is not asymptotic and does not suffer from the sequential 
dependence issue explored in HMoi Let the deformation be y(x) = Fx, x E £1. The energy equals 
the exact sum of the elastic energy J Q W(F)dx plus the surface energy j m ^ (F,n)dx, plus the 
corner energy J2iLi r (-^' n «' n i-i)> summed over the N vertices of fl. The surface energy density is 
explicitly obtained: 

70(^*0 = -t E ^(K^i-iMI^I), (1.2) 

4 / — ' n 

weL\{o} 1 1 

where n is a normal to dVl whose components on each facet are the Miller indices (irreducible 
integers) of the corresponding lattice plane, and <p is the interatomic potential. The corner energy 
t(F, rii, rii-i) is also explicit but more complicated; apart from F, it depends on the two unit normals 
of the facets meeting at the ith vertex. 

For an infinite range potential this result retains only asymptotic validity for a lattice polygon kVt 
as k — > 00; the three energies just mentioned are the first three terms of the asymptotic expansion 
of the energy for large k (Proposition [33]). The order of the remainder depends on the decay rate of 
the potential (p. 

Because of its construction based on lattice polygons, the surface energy density ( 11.2b is only 
defined for "rational" directions of the surface normal (n = (ui, v 2 ) E S 1 is called rational if u 2 / V\ 
is a rational number or V\ = 0, irrational otherwise). It is natural to ask how (11.2b can be extended 
to irrational normals. This is important because when f2 is strictly convex and OVl is smooth for 
example, the normal to d£l is irrational almost everywhere on d£l. We study this issue in Section HI 
The main idea here is to approximate such a strictly convex, smooth f2 by a suitable lattice polygon. 
There is a natural way to do this: we consider the convex hull of all lattice points contained in 
f2. This allows us to use number-theoretic results on the asymptotic properties of such sets due to 
Barany and Larman UBIi and also Huxley HHul ; see also the survey [IKMJ. Perhaps surprisingly, the 
surface energy density we obtain for smooth strictly convex domains (Proposition 14. lb is different 
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from (11.21) . It takes the form 



1 



\w-n\<p(\Fw\), 
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(1.3) 



w£L\{0} 



where n is now the unit normal to dQ and can take on irrational values. This allows us to extend 
the definition of the surface energy density (11.21) to all values of the normal. We consider more 
general domains that are convex but whose boundary comprises flat facets as well as curves with 
positive curvature. In Proposition 03] we calculate the surface energy of such domains. It equals the 
boundary integral of an extended surface energy density function, 



with 7<> from (11.21 ) and 7 D from (11.31) . The dependence of the surface energy density on the surface 
normal is thus quite intricate. One can show that ^(F, •) : S 1 — > R is continuous at irrational n, 
discontinuous at rational n, and almost nowhere differentiable (Proposition 14.41 ) . Because of the 
above pathological properties, the surface energy density does not satisfy the usual hypotheses of 
the Wulff theorem; see e.g. [Fo|. 

This difficulty is dealt with in Section [5J One observes that given a collection of lattice points 
within a convex domain, there is some freedom in choosing an alternative convex domain containing 
precisely the same lattice points. The choice of this domain plays a subtle but important role in the 
description of surface energy. Utilizing this, we are able to write the surface energy of the crystal in 
the form 



where Q' is a suitably conrtructed domain containing the same lattice points as fl The surface 
energy density 7 D , defined in (1 1 .3ft , is actually Lipschitz continuous in the unit normal. 

For finite range potentials, Herring HHeH and Parry llPal give formulae for the surface energy of 
a planar surface that bear some similarity with 7 Q of (11.31) . 

A more realistic approach to surface energy would allow for local "relaxation" of atomic po- 
sitions from the macroscopic deformation near the boundary. Such deviations would perhaps be 
determined by minimization of the atomistic energy. This seems to be a formidable problem in the 
setting considered here (more than one dimension, general boundary geometry, arbitrary interac- 
tion range, nonconvex potentials); there are very few explicit calculations of such a relaxed surface 
energy, even in substantially simplified settings. One of the first results in this direction is due to 
Braides and Cicalese HBCII ; they obtain the relaxed surface energy of a one-dimensional atomic 




lo(F,n) 
lo{F,n) 



n rational (n/|n| = n) 
n irrrational, 



(1.4) 
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chain using T-convergence. The result is not explicit and seems difficult to compare quantitatively 
with the explicit "constrained" energy of Mora-Corral [MoJ. In a two-dimensional setting, Theil 
HTM calculates the relaxed surface energy of a crystal with quadratic short range potentials; the 
result is in the form of a perturbation of the constrained surface energy explored here. 

In order to obtain some quantitative information on the difference between the relaxed and con- 
strained surface energies, numerical optimization of the atomistic energy was recently performed for 
a completely unconstrained, Lennard- Jones two-dimensional crystal [RoJ. Atomic positions were 
allowed to relax from initial positions forming a lattice triange or hexagon with low Miller-index 
boundary. The constrained energy was obtained by minimizing over the deformation gradient ma- 
trix of a homogeneous deformation that the atoms are constrained to follow. It was found that the 
difference between the relaxed and constrained surface energies is typically less than three percent 
(after the appropriate scaling and bulk energy is accounted for). This suggests that there are some 
situations where the constrained surface energy is close to the unconstrained one. In analogous one- 
dimentional computations, the results agree qualitatively with the conclusions of HBCH : at the same 
time, the difference between relaxed and constrained surface energy is less than one percent. 

Many of the results presented here, in particular expressions (11.21) through (11.41) for the surface 
enrgy density, are valid for three-dimensional crystals as well URol . The three-dimensional version 
of Pick's Theorem that these results rely upon is not completely settled at this point HBPl fBRI. 

2 The Bond Counting Approach 

For subsets P, Q of W 1 , define the Minkowski sum P (B Q = {p + q ■ p E P, q E Q} and write 
p + Q = {p} © Q. The lattice is L = Z 2 unless otherwise noted. This is not a restriction; the results 
can be adapted to any Bravais Lattice L* by incorporating the linear mapping from L onto L* into 
the deformation. For x = (a, (3) £ Z 2 let 

gcd(x) = gcd(|a|, \/3\), x= 1 x, x ± = ((3,-a). 

gca{x) 

We assume that the reference region fi C K 2 is a convex body, or a compact convex set with 
nonempty interior. One important case we will consider is when C IR 2 is a convex lattice polygon. 
In particular, = convj^i, . . . , vn}, the convex hull of its iV vertices Vi E L, i E {1, . . . , N}, 
which are lattice points. The boundary dVl consists of N facets Si = conv{t>j, v i+ i}, where v^+i = 
v\ and Sn+i = S\. Letting m.j = v i+ i — V{, fhi = rxiij gcd(mj), the Miller normal to Si is 
fii = rhf, so that gcd(nj) = 1. Fix w E L, let x E L and define b = b(x,w) = {z E M 2 : z = 



Surface Energy 7 

x + tw, < t < 1} = conv{x, x + w} as the bond starting at x with bond vector w . The set 

B W (Q) = {b:b = b(x,w),x E Q H L, x + w E Q H L} (2.1) 
is the set of all w-bonds of f2 (bonds with bond vector w). We will use the abbreviation 

60 = 6(0, w). 

The energy of the homogeneous deformation = Fx can be written as 

E{n, y }= 1 - J2 E vd^D- < 2 - 2) 

x'GQnL w£L\{0} 
x+wGQnL 

The factor of 1/2 occurs since b(x, w) = b(x + w, —w) and the potential (p is even in w. Interchang- 
ing the order of summation above we obtain 

m,y} = l E E v(i^H) = ^ E E v(i^i) = i E E L 

«>eL\{o} xef^ toei\{o} beB w (n) weL\{o} beB w (n) 

x+w€QnL 

Evidently, in order to determine the energy, it suffices to calculate, for each w E L, the w-bond 
number of Tl, i.e., N W (Q) = #B w (Vt); see (12.11) : 

E{tt,y}= 1 - f(\Fw\)N w (n). (2.3) 

w&L\{0} 

Clearly the number of w-bonds "starting" in equals the number of lattice points of fh 

#{6 = b(x, w) : x e n n L} = #(n n L). 

Some of these bonds are not contained in 5 TO (0): 

N w (n) = #(anL) -#r w (fi), T m (0) = {6 = &(i 1? D):iGfinL,i + «j^ni}. (2.4) 

Recall that 5 = <9fi consists of iV facets Si, i = 1, . . . , N, each with unit normal rii, outward with 
respect to f2. Let 

J(w) = {i e 7L : 1 < i < N, ^ ■ w > 0} , = |J S u (2.5) 

ieJ(u>) 

so that is the part of <9f2 through which u> points outwards. Denote by the set of all 

w-bonds that intersect and terminate outside f2. 

T^(fi) = {6:6 = 6(x, iu) E B W (L), 6 n S+ ^ 0, x + w £ O}. (2.6) 



8 



Rosakis 



Some of these bonds "straddle" Q, that is, have both endpoints outside Q but intersect dQ; specifi- 
cally, 

= {b(x, w) e Tl{Q) : x fl, x + w £ O}. (2.7) 
Then obviously in view of (12.41) . 

T w (n)=Tl(Q)\Tl(Q). 

As a result, 

N W (Q) = n L) — #T+(^) + #T*(Q). (2.8) 

Roughly speaking, the number of w-bonds in equals the number of lattice points in it, minus the 
number of bonds that traverse the boundary at least once, plus the number of bonds that traverse the 
boundary twice. The reason for the splitting (12.81 ) is that each term can be evaluated using results 
from geometric number theory. 

The number of lattice points in n,#(flni), is addressed by Pick's Theorem, IIPTllRell , a variant 
of which is the following 

Lemma 2.1. Let Q be a simple closed lattice polygon with facets Si and outward Miller normal 
n = fti on Then 

1 - 191 

#(finL) = |fi| + -Vp[ + l. (2.9) 

1=1 

Equivalently, letting 9i be the (dihedral) angle between normals of facets meeting at the ith vertex, 

#(nr\L)= f idx+ f -Lds + T^. (2.10) 

Jn Jan 2|n| 2vr 

Proof. Pick's Theorem [|PillRell states that 

|q| = #(n n L) + ^#(<9ft n L) - 1 = #(Q n L) - ^#(<9ft ni)-i (2.1 1) 

(since is closed). If two neighboring lattice points in a facet 5, differ by m s G L (with relatively 
prime components), then n L) = \Si\/\rhi\ + 1 while #(<9fi n L) = Y*Li[#( s i n L) — 1] = 
|5i|/|mj| since each Si contains both its endpoints. Now the Miller normal n,j = fhf, so that 
\n%\ = \rhi\ and (12.91 ) follows. Also, (|2.10l) is a trivial consequence of (12.91) , given that the sum in 
(I2TT01) equals 1. □ 

Eq. (12.101 ) has an interesting interpretation. It exactly equates a discrete quantity (number of 
atoms in f2) with a continuum expression: the "volume" integral of a bulk density, plus the "surface" 
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integral of a surface density, plus contributions of corners. We will show in the sequel that both the 
w-bond number N w (Cl) and the energy admit analogous representations. 

The first term in (|2.8I) is given by (12.91) . Turning to the second term, let Pj (w) be the parallelogram 
60 © Si with two parallel sides Si and w + Si if w ■ rii > 0, Pi{w) = otherwise. Then it is easy to 
see that b(x, w) G T^(O) if and only if x + w G Pt(w) \ Si for some i G J(w); see (12.51) . Thus 

2j(n) = {6(ar,u;):x + «;GP(«;)ni}, P(w)= |J P^w) \ Si = (b © 3+) \ S+. (2.12) 

ieJ(w) 

It follows that 

#2l(n) = #(p(«;)nL). (2.13) 

In general, P(w) is not convex. However, if one defines 

n w = b ®n= \J(tw + n), (2.14) 

te[o,i] 

then Cl w is a convex lattice polygon, being the Minkowski sum of two such sets. In fact, 

£l w = conv{fi, w + fi}. (2.15) 
Also P(w) = \ f2, while VI C £l w . This and (12.131 ) imply 

#Tji(0) = n L) - n L). (2.16) 

The right hand side can be evaluated using Lemma 12.11 for each term. Note that dQ w comprises 
d£l \ S+, w + S+ and two w-bonds joining these two pieces. The result is 

N 

#Tt(fi)= \S l \w-n i + \b\/\w\ = J2\S l \(w-n l )+g C d(w), (2.17) 

ieJ(w) i=l 

where (x) = (x + \x\)/2 for s G R and n = n.j on Si is the unit outward normal to <9fi. Here 

= gcd(w). 

It remains to evaluate T|(f2). If a bond 6 = w) terminates in w + f2, or x + w G w + f2, 
then x G fl This together with (12.71) and (|2.12l) immediately shows that b G T^(O) if and only if 

x + w G P(u>) \ w + VI. Since P(u>) = Q W \Q, 

= #(QH n L), qw = fi w \ (n u (™ + n)). (2.18) 

We will show next that for |iu| small enough compared to the facets of Q, Q(w) consists of one 
or two triangles, each having a vertex at one of the two ends of the simple polygonal line S^. For 
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example, if f2 = [0, 3] 2 and w = (1, 1), Q(w) consists of the triangle with vertices (0, 3), (1, 4) and 
(1,3) and its image under reflection about the (1, l)-axis. Any b 6 Tjj(fi) intersects two different 
facets of dVl by (|277T) . Let 

5 = 5(Sl)= min dist(^, Sj) (2.19) 

l<i,j<N 
Vi(ZSj 

where Vi E Z 2 are the vertices of f2. The shortest line segment with endpoints on non-adjacent 
facets has length 5. If |iu| < 5, b 6 T$(Q) necessarily intersects two adjacent facets, say Si and Si-i 
meeting at some vertex with outward normals rii, n^i (where n = n^). Since both endpoints of 
b are outside fi, w ■ rii and w ■ n^i must have opposite signs. Then in case w ■ rii > and w ■ n^i < 0, 
x + w is in the triangle with vertices Vi, Vi + w and the intersection of S{ and w + which is 
therefore is part of Q(w). If the reverse inequality holds, the triangle with vertices Vi, Vi + w and the 
intersection of w + Si and S^i is part of Q(w). Regarding lattice point count, both cases reduce to 
the triangle with base b and sides normal to n, and r^_i: 

r(w,nj,nj_i) = conv{0,w,g}, g-n; = 0, (g- w) -nj_i = 0, (u>-7ij)(«;-nj_i) < 0. (2.20) 

In addition, the relative interior of the base b(vi, w) of the triangle with endpoints t> , + w is also 
part of Q(w) and contains gcd(u>) — 1 lattice points. Consequently, if \w\ < 5, 

#Tl(Q)= [gcdH-l + #T(^,n i ,n i _ 1 )]. (2.21) 

l<i<N 
(wni)(w-m-i)<0 

Unfortunately, T(w, rii, rij_i) is not a lattice polygon in general, since q need not have integer coor- 
dinates and Lemma [27T1 does not apply. Instead, we count the lattice points inside the triangle more 
directly: 

Lemma 2.2. Suppose (w ■ rii){w ■ n^-i) < and let T = T(w, rii, n i-i) C R 2 be the triangle of 
(12.201) . Let ti 6 Z 2 be such that {u, w} is a lattice basis for Z 2 . Then 

#(T fli) =JV r (tw,n i ,n i _i), 

where for tcfZ 2 and unit n, m G R 2 with (it; ■ n){w ■ m) < 0, 

777/ • 

Nt(w, n, m) = s (q ■ w x , q ■ u x , gcd(w)) , g = -n' L (2.22) 

m ■ n 1 - 

and s : R x R x Z R is given by 

s(a,P,k) = (l-\\a\])(k-l)+ l ^T (\^±j\- f-j) (2.23) 
with s(0, /3, Jfe) = 0. 
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Proof. Let n — rii,m — 7ij_i. Since in (|2.20l) q ■ n = 0, q = Xn 1 - for some A G 1. Then solving 
(q — w) ■ m = for A gives q as in the second of (12.221) . Let w = w/ gcd(w) = (wi,w 2 ) and 
suppose u — (ui, u 2 ) G Z 2 solves u ■ w 1 - = 1, or w 2 Ui — w 2 u 2 = 1. This is solvable by Bezout's 
Lemma since gcd(wi,w 2 ) = 1. Then the matrix A = co\(u,w) has unit determinant u ■ w 1 = 1 
and integer entries, hence so does A^ 1 = row(w ± , u x ). As a result {u, w} is a lattice basis for Z 2 , 
while the linear transformation with matrix A^ 1 is lattice invariant . Now T = A~ l T has vertices 
0, (0, k) G Z 2 and p = (a, (3), where 

fc = gcd(tw), (a,(3) = (q-w ± ,q-u ± ); (2.24) 

in general p is not a lattice point. Suppose for the moment that a > 0. Then 

/3 , 0-k 



f'nl = < (xi,x 2 ) G Z 2 : < X\ < a, — Xi < x 2 < k + 
{ a 



a 

For x G R let [x\ ' be the greatest integer strictly less than x and |~af| ' the least integer strictly greater 
than x. Then the number of lattice points on a segment {(x 1 , x 2 ) : x± = j, p < x 2 < u}, where 
j G Z and /i<i/GR equals [^J' — [//]' + 1. Hence, 

H' 



#(f nL) = £( + — J '- + 



i=i 

Since [xj' = [x] — 1 and |V|' = [x\ + 1, the above reduces to s(a, (3, k) in (12.231) . It then follows 
from (12.241) and (|2.22l) that #(T' fl L) = Nt(w, n, m). The linear transformation with matrix A is 
lattice invariant and thus #(AT' n L) = #(f ' n L) llEFl . while AT' = T. In case a < 0, reflect T 
by replacing a by |a|. If a — then T" = 0. □ 

This together with (12.211) gives 

#T*(ft)= ^ [gcdH-l + iVrK^ni-i)]. (2.25) 

l<i<Af 
(w-m)(w-ni-i)<0 

To obtain an expression for the u»-bond number of f2, merely substitute (12.91) , (12.171) and (12.251) into 
(12.81) and rearrange. Observe that for a given bond vector w, N W (Q) is completely determined by 
the area the lengths \Si\ of the facets, and their orientations through the Miller normals nf. 

Lemma 2.3. Suppose \w\ < 5, cf. (12. 191 ). Then the w-bond number of is given by 



1 N ISA 

(ft) = \n\ + -^ l -^(l-2(w-ni))+l-gcd(w)+ Yl [gcd(w)-l+N T (w,n i ,n i ^ 1 )}. 

i=l ' l<i<N 

(2.26) 
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The above can readily be written in a form similar to (12.101) : 



N W (Q)= / ldx + / g(w,n)ds + y^ j h( 
Jn Jan i=l 



w,ni,ni-i), w e L, 



as a bulk integral, plus a "surface" integral, plus corner contributions, for suitable normal-dependent 
densities g and h. 

The present approach of counting bonds has certain similarities with the bond density lemma of 
Shapeev (Shi. 



3 Surface Energy of Lattice Polygons 

We are now in a position to compute the energy. Consider first a finite-range potential that only 
involves bonds within a bounded set. Let the bond range R C L\ {0} be symmetric, so that w E 
R =>■ — w E R. Allow the interatomic potential (p w (-) to depend explicitly on w, require 

ip w {-) = ip- w {-) \/w E L, (p w (-)=0VweL\R, (3.1) 

and define the energy of the homogeneous deformation y(x) = Fx, x E f2, 

xeQflL w£R 

where ip w : (0, oo) — >■ R is not restricted to be regular in any way. 

Proposition 3.1. For F E M 2x2 , fh E Z 2 and unit n,m el 2 define the stored energy function 

W{F) = \Y,V w {\Fw\) ) (3.3) 



2 

weR 



the surface energy density function 



7o(F, m) = ~Vi (|w • m| - 1) <^(|Fu>|) (3.4) 
4 ^— ' m 

and the vertex energy function 

t(F, n, m) = \ I H n<m {w) - ^( n , m ) (gcd(w) - 1) + H n , m (w)N T (w, n, m) j <^(|Fw|), 



(3.5) 

where the sector step function 



H n , m (w) 



1 if (if • n)(iu • m) < 0, 
if (w ■ n)(w ■ m) > 0, 
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and 9(n, m) is the angle between n and m, while N T is defined in Lemma \2.2\ Suppose the bond 
range R is bounded with max weR \ w\ < 8, cf. (12. 19b . Let n = fti on Then the following expres- 
sion is exact: 

E{VL,y}= [ W(F)dx+ [ ^(F, n)ds + YV(F, n u m_ x ). (3.6) 
Jn Jon i=1 

Proof. As in the argument leading to (12.31) . one can write (|3.2I) as 

E{Q,y] = ^^2N v) (Q)( Pvi {\Fw\). 

By the hypothesis on R, Lemma [231 holds for all w G R. Multiply (12.261 ) by ip w (\Fw\) and sum the 
result over w G R. Interchange the order of summations, noting that 

^2(w ■ ni)<f w (\Fw\) = - ^2 \ w ■ ni\(f w (\Fw\) 

uiG-R w£R 

by the symmetry of R and the first of (13.11) . also that the sum of the (dihedral) angles between 
normals of facets meeting at vertices J2 i=1 0(rii, n^i) = 1, and finally that summation over w in 
the sector of R where (w • n)(w • m) < can be replaced by summation over R provided the 
summand is multiplied by H n ^ m (w). □ 

Next we consider infinite-range potentials, where R = L \ {0}. We seek the energy of the kth 
dilation kVt of the domain, k G Z + . Here we have no choice but to let k k be an integer; otherwise 
kVl is not a lattice polygon in general. The following will be useful. 

Lemma 3.2. Let M be a positive integer. For p > sufficiently large and p > 0, 

\w\- (M+p) < c P - p 

w£Z M , \w\>p 

(where C > is independent of p and p). 

Proof Let x : ~R M — > 7h M be the map x(J2fLi x i e i) — Yl!^=i\. x i\ e ^ x -i e ^> w i m e « standard basis 
vectors for IR M . Thus \x — x(x)\ < D, the unit cell diameter. Write x(x) = x + (x(x) — x) and 
invoke the triangle inequality to conclude \x\ — D < \x(x)\ < \x\ + D. Then also \x(x)\ > p 
implies \x\ > p - D, while |x(x)|~ (3+p) < ||x| - D|" (3+p) . Thus A p = {x G R M : \x(x)\ > p} C 
R M \B P „ D (0). 




r oo roo 

<C (ar)-( M+p) r M - l dr = C r - {1+p) dr = Cp- p , 

J ap J p 
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where a G (0, 1) is such that ap = p — D, so that a G (1/2, 1) for p > 2D and C takes possibly 
different values on different sides of an equation. □ 

For convenience we suppose that the interatomic potential ip w (-) = ip(-) (does not explicitly 
depend on w), although this is not essential. 

Proposition 3.3. Suppose the interatomic potential ip : (0, oo) —> K satisfies the following: for each 
r > and for some constants C = C(r ) and d > 2, 

\ip{r)\ < Cr- {2+d) forr G [r„,oo). (3.7) 

Let the bond range R = L \ {0} in (13.21 ) and in the definitions (13.31 ) ofW, (13.41) 0/70 and (13.51) of 
r. Then as k — > 00, k G Z + , 

E{ktt,y} = k 2 [ W{F)dx + k [ 7,(F,n)rfs + Vr(F,fi ! ,)i ! _i) + 0^"' i ). (3.8) 
Jn Jan i=l 

Proof. Note that 8{kVL) = k5(Q) = kb in (12.191) . so that Lemma[23]fbr kVL holds provided 

weR k = (L\ {o}) n 5 M (o) = (L n s M (o)) \ {o}. (3.9) 

Split the energy as follows: 

E{kn,y} = ^Y1 N w (kn)p(\Fw\) + i ^ iV tt (A;n)^(|F«;|). (3.10) 

Now it is clear that for any w E L and A; G Z + , 

< A^(fcQ) < #(£;fi n L) < Ck 2 

for some constant C > 0, since all bonds within A;fi start in kVt and by Lemma I2TT1 applied to kVL 
(the dominant term in (12.91) would be \kVt\ = k 2 Vt). This provides a bound for the second term in 

dug) : 



w&L\R k 



<C£; 2 ^ ^(1^1)1 <C£; 2 ^ |cH~ (2+d) < C£; 2 - d , 

w£L\R k w&L 2 , \w\>k8 

(3.11) 

where we invoked (13.71 ). a > is such that > a\z\ for all 2 G M 2 and we used Lemma I3T21 with 
p = k5 and p = d; C is a generic constant with different values on either side of an (in)equality. 
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The first term in (13.101) is covered by Proposition 13.11 applied to kCl, since w E R k means 
\w\ < kS = 5(kQ). Noting that \k£l\ = k 2 \tt\, \kSi\ = k\Si\, Proposition ^ . 1 1 implies 

. N N 

- ^(^^(iF^^fc^l^^ + fc^l^lTfe^nO + ^r^F,^,^.!), (3.12) 

wGRk 1=1 i=l 

where W k , ^ k and r k are given by (13.31) . (13.41) and (13.51) with in place of R; see (13 .9b ■ Recalling 
that W, 7 and r are defined by the same equations with R = L \ {0}, using Lemma [3T2l with 
M = 2, we may estimate (omitting arguments) 

|W- < CAT d , |7-7 fc | < Ck l -\ \r - r k \ < Ck 2 ~ d . (3.13) 

We only demonstrate the third of these, the others being easier. Recall that in (13.51) . Nt is the 
number of lattice points in the interior of a certain triangle T whose area is bounded above by 
C|w| 2 , cf. Lemma [2T2l By Pick's Theorem (12.1 II) (applied to the lattice parallelogram of smallest 
area A containing T, and having the same base) the area A exceeds Nt hence Nt < C\w\ 2 . Also 
gcd(w) < |u>|, \H nym \ < 1, hence we have from (13.51) . 

\r-r k \< C\w\ 2 \ip{\Fw\)\< C M~ d < Ck 2 ~ d 

w)£Z 2 . \w\>kS w^X 2 , \w\>kS 

proceeding as in (13.111) . By (13.131) . replacing W, 7 and r by W k , ^ k and r k in ( 13.121) produces an 
error of 0{k 2 ~ d ). Combine this with (|37TTT> and (13T01) to obtain CT) . □ 

4 Surface Energy for More General Boundaries 

We examine the surface energy density function 7 in (13.41 ) more closely, paying attention to its 
dependence on the surface normal. Due to its construction, 7<>(F, •) : M — > M. is defined only for 
"rational directions", that is, on the set of Miller normals 

M = {n:n= (is u v 2 ) G Z 2 , gcd(z/ 1; u 2 ) = 1}. (4.1) 

Using (13.31 ), we rewrite 7 in (13.41 ) as 

^(F,n) = --^\w -n\<p(\Fw\) + 2^-W(F), n = n/\n\, neM. (4.2) 

The first term (involving the sum) reduces to a function of the unit normal n, and trivially admits a 
unique continuous extension onto the whole of the unit circle S 1 . There is no such extension for the 
second term. Define the rational and irrational direction sets as 

S\ = {n:ne S\ n = n/\nl ne M}, S] = S 1 \S l R , (4.3) 
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respectively, where M is defined in (14.11) . Thus a vector is rational (irrational) if the tangent of the 
angle it makes with the usual basis vectors is rational (irrational). Since facets of lattice polygons 
have rational normals, the surface energy density 7<> is defined only for such directions. Note that 
for each n E Sr there is a unique n = n(n) E M with n/\n\ = n. The question arises as to how one 
can extend the definition of %(F, n) = y (F, \n(n)\n), n E S^, to the whole of S 1 . This is related 
to another question: what is the surface energy when d£l is smooth, for example d£l = S lc ! It turns 
out that this question can be answered, at least partially, using the present approach. The basic idea 
is that even if dVl is not polygonal, but smooth, the convex hull of all lattice points inside is a 
convex lattice polygon. 

Proposition 4.1. Let Q C R 2 be strictly convex and dVt be C 2 with positive curvature. Suppose ip is 
as in Proposition ^. 3\ but with d > 3. Define the reduced surface energy density y Q : M 2x2 x S 1 — > R 
by 

j (F,m) = -- \w-m\(p(\Fw\), F E M 2x2 , mES 1 . (4.4) 

w£L\{0} 

Then for any sequence r = r k — oo as k — oo (r k E R + , k E Z + ), 

E{rQ,y} = r 2 [ W(F)dx + r [ j {F,n)ds + 0(r 2/3 ), (4.5) 
Jn Jan 

where n : dQ — > S 1 is the unit outward normal to dQ. 

Proof. For each r > 0, let f2 r = conv(rfi fl L). Then Vt r C rVt is a convex lattice polygon, while 
rflnl = fl r nl. Hence, in view of (12.21 ), 

E{rn,y} = E{n r ,y}, (4.6) 

where y(x) = Fx for x E rQ. The calculation of E{VL r , y} proceeds as above with one exception. 
For any w E L, the condition |iu| < 5(Cl r ) (see (12.191 )) may be violated for large enough r, since 
facets may become as small as the shortest lattice bonds. This affects N w (Cl r ), but only the part 
regarding #T^(f2 r ) — see (12.81 ), (12.181) — which we merely need to estimate. Given any convex body 
D C R 2 , let 

Q W {D) = (bo © D) \ [D U (w + D)], Q W {D) = Q W {D) n L, 

where b = conv{0,w}. Then by (TTT81) and (I27T5T) . #T|(fi r ) = #Q(fi r )- Since tt r c rQ, and 
rfi \ £l r contains no lattice points, (12.181 ) and (12.151) imply 



Q w (tt r ) C Q w (rft). 



(4.7) 
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Let q, q' C d(rQ) be the two points of d(rQ) where the tangent vector is w, and B rp C rf2 be a 
disk with <92? rp tangent to <9(rf2) at g, where p is the smallest radius of curvature of dVt. Also let 
B' C rO be a similar disk tangent to d(rQ) at q'. Then for r large enough it is easy to see that 

Q w (rty C Q w (B rp ) U Q w (B' rp ). (4.8) 

The connected component of Q w (B rp ) containing q is contained inside an isosceles triangle with 
base a w-bond (with length \w\), and height the distance from the base middle to the intersection of 
the two circles dB rp and w + dB rp ; these are tangent to the base at its endpoints. The triangle height 
is thus bounded by C(r)\w\, where C(r) approaches zero for large r. A crude but sufficient upper 
bound of the lattice point count of this set, hence also of the right hand of (14.81) . is C|w| 2 , with C 
independent of r. In view of of (|4.7I) . #T|(fi r ) is also bounded by C| w| 2 . This estimate replaces the 
sum over vertices (second sum) in (12.261) . Since J2weL\{o} \w\ p (f(\Fw\) are absolutely convergent 
for p — 0, 1, 2 as one infers from Lemma l3?2l it follows that 

E{n r ,y}=\n r \W(F)+ [ 7o (F,n)ds + 0(l) 

\tt r \+ [ -^rds W(F)+ [ lo {F,n)ds + 0{l) 

=# (rfi fi L) W(F) + / 7o(F,n)tfe + 0(l). (4.9) 

Jdn r 

Here we have used (14.21) and (14.41) . then (12.101) . in which the last term (sum) equals 1, together with 
the fact rfi n L = (l r fl L. We turn to J m J (F, n)ds. Recalling (14.41) . a typical term involves 

\w n\ds = 2\Pro) w ^dVt r \\w\ } (4.10) 

\Pro) w ±dVl r \ being the length of the projection of <9f2 r onto a line perpendicular to w. This follows 
after splitting dVL r into two pieces, over which w ■ n is > and < 0, and using the Divergence 
Theorem on each. Next, we show that 

< |Proj w xd(rfi)| - |Proj w x9fi r | < C\w\ 2 , (AM) 

where the constant C is independent of r > 1 and w. There are lattice points z~ and z + G dVL r D 
L, such that f2 r lies entirely between lattice lines /~, / + with normal w 1 - and containing z + , 
respectively. Consider the part of <9(rf2) that lies outside the strip bounded by / + and l~ . It consists 
of two disjoint arcs, one to the "right" of Z + and the other to the "left" of l~ . The length of the 
projections of these two arcs onto the w 1 - axis equals the difference in (14.1 II) . Let c + be the arc to 



18 



Rosakis 



the right of l + (with endpoints in / + ). Let s be the region bounded by c + and l + . The only lattice 
points it contains are in Z + . This is true since rQ \ £l r is free of lattice points. By the strict convexity 
of rf2, there is a unique q E c + where the normal to c + is w 1 . Consider the osculating circle of c + 
at q. Let s' C s be the portion of the osculating disc contained in s; it is a circular segment whose 
height (in the direction w 1 -) equals the thickness of s (the length of its projection onto a line along 
w 1 ). The radius of the circle is rp for some p > 0. There are two possibilities. Either s' lies between 
l + and the next lattice line V with normal w to the right of l + , or it extends beyond I' to the right. 
In the first case the height of the segment s' is l/|w|, the distance between adjacent lattice lines 
with normal w^. In the second case, let s" be the portion of s' to the right of /'. Then s" is also a 
circular segment and free of lattice points. Suppose its chord length is c and height is h. Since the 
radius of the circular arc is rp we have h 2 — 2rph + c 2 /4 = 0. Solving this for h/{rp) and using 



the inequality 1 — y/1 — x < x for < x < 1 yields h < c 2 / (4rp). Now since the circular segment 
s" is free of lattice points and its chord is in I', the chord length c < \w\ < \w\ (since the distance 
between adjacent lattice points in V is |u>|.) Hence h < \w\ 2 /(Arp). The total height of the larger 
circular segment s' is h+ l/\w\ which is thus bounded by C|w| 2 for r > 1. The thickness of s in 
the direction normal to w is the same as this height. This shows (14.1 II) . 
Combining (14.111) with (14.101) shows 

<C\w\ 3 . (4.12) 



\w ■ n\ds — \w • n\ds 

8(rU) Jdn r 

In view of (14.41) and since the sum J2 w &l\{o} \w\ 3 (p(\Fw\) converges absolutely by hypothesis, one 
deduces 

/ lo (F,n)ds= [ lo {F,n)ds + 0{l)=r [ 7o (F, n)ds + O(l). (4.13) 

Jdn r Jd(rQ.) Jdn 

Our hypotheses on dQ ensure that #(rf2 D L) = r 2 \Q\ + 0(r 2/3 ), e.g. llHulUKMi This together 
with (I4T3T) in (EU) and g3 proves (14.51) . □ 

According to Proposition 14. 11 when <9f2 is smooth and strictly convex, so that the normal vector 
is irrational almost everywhere on <9f2, the surface energy density is given by (14.41 ); in contrast, for 
lattice polygons (with rational normal a.e. on dQ), the surface energy density is given by ( 14.21) . This 
suggests that we combine the two expressions in defining a surface energy density for all values of 
the unit normal. That will allow us to treat the general case when £1 is a (not necessarily strictly) 
convex body. We do place some restrictions on d£l: flat parts of dQ must be lattice segments (with 
rational normals). Corners have to be lattice points. 

We will need the following auxiliary result: 
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Lemma 4.2. Let D C M 2 be a strictly convex body and <9-D be C 2 with positive curvature. For 
r > define the convex lattice polygon D r = conv(rD n L) with Miller normal n : <9-D r — > M. 
Then as r — > oo, 

/ _L ds = 0(r 2 / 3 ). 

Proo/: By Pick's Theorem (Lemma Q>, and since #p r nl) = #(rL> n L), 

[ J-ds = n L) - \D r \ - 1 = #(r£> n L) - \tD\ + IrDl - \D r \ - 1. 

■W 2|n| 

Now #(r£> n L) - |r£>| = o(r 2 / 3 ) by JHul. In view of Theorem 4 and Remark 2 of [BLJ, and since 

D r C tD, 

< |r£>| - \D r \ < Cr 2/3 (4.14) 
for some constant C. The result follows. □ 

We now state the main result of this section: 

Proposition 4.3. Assume that Q is a convex body with dQ Lipschitz, and that there is a finite set 
of lattice points {v%, . . . ,v^} C dQ R L, that partitions dQ into N curves Si, dQ = (J i=1 each 
with endpoints Vi and v i+ \ (v^+i = Vi), such that Si fl S i+ i = v i+1 , Si is a C 2 curve and one of the 
following two alternatives holds: 

(i) For i G J c C {1, . . . , iV}, S{ C Tj, where Ti is a simple closed C 2 curve with positive curvature, 
or 

(ii) Fori G J/ = {1, . . . , iV} \ J c , Si is a straight segment. 

Suppose if is as in Proposition \3.3\ but with d > 3. Define the extended surface energy density 
j(F, •) : S 1 -»■ R as follows: 

1 1 - 

■- ^2 \ w ■ n \<f(\ Fw \) + ^T~i W ( F ), n e S l R (n G M, fi/\n\ = n), 

w£L\{0} 



l(F,n) = < 



(4.15) 

~I \ w ' n \ ^(l-^H) = lo{F,n), neSj, 



4 

u;£L\{0} 



with 70 defined in (14.41) and Sr, Sj defined in (14.31) . Then as k — > oo, k G Z+, 

E{k c l, y} = k 2 [ W(F)dx + k [ ^(F,n)ds + 0(k 2/3 ), (4.16) 
in Jan 

where n : <9f2 — )■ S* 1 ?s unit outward normal to dVt. 
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Proof. We now choose r = k E Z + and let Q k = conv(kQnL). The part of the proof of Proposition 
14. II prior to (14.91 ) is easily adapted to the present setting, so that once again, as k -» oo, with 7o as 
infill). 

£{fi fe ,y} = |^|^(F) + / lo (F,n)ds + 0(l). (4.17) 



Let <9f2f be the union of those Si that are straight segments and dVl c the union of the Si with positive 
curvature, so that dVl = dVtf fl dVt c . By hypothesis, for k E Z + we have kvi E d(k£l) fl L, hence 
also kvi E d(Q k ) n L. Then fcdfi/ C dQ k = d(Q k ). Let <9fi£ = dQ k \ kdQ f . Then 

E{n kl y} = \Q k \W(F)+ [ ^(F,n)da + / 7<> (F,n)ds + 0(1). 

Our hypotheses regarding dVt c , specifically alternative (i), ensure that n E S] a.e. on kdVL c , while 
(ii) implies that n E S] a.e. on kdVLf. Using (14.151) , rewrite the above as 



E{n k ,y}= \m\W(F)+ [ lo {F,n)ds+ f j(F,n)ds + R(k), (4.18) 

Jkdn c Jkdflf 

where 

R{k) -- 



\ft k \ — \kQ\ + I ——ds 



W{F)+ [ lo {F,n)ds- [ 7o (F,n)ds+0(l). (4.19) 
Jan? Jkdn c 



It remains to show that R(k) = 0(k 2 ^ 3 ) as k — > oo, k E Z+. Let i G J c , so that 5j satisfies 
alternative (i) in the statement of Proposition 14.31 Let S k be the portion of dVt c k between kvt and 
kvi + i, i.e., terminating at these two points and containing no other kvj. Let the strictly convex 
body D l be such that dD l = Tj. Let G\ be the bounded region whose boundary is kSi U S k ; 
this is well defined since both curves terminate at kvi and kv i+ \. Then G\ C kD l \ D l k , where 
Dl = conv(kDi fl L), and 

\kn \n k \ = j2 \g{\ < \ kr>i \ D i\ < Ck2/3 ( 4 - 2 °) 

in view of (14.141) applied to D l for r = k E Z + . 
Next, note that Sj(. C <9.DjL As a result, 

< / -^cfc = V / -t-ds < V / -^ds < CA; 2/3 (4.21) 

/flo c 2 n z — ' /<?; 2 n ^— ' lani 2 n 
• /c ' w fe 1 1 ie,j c Jb k 1 1 i£j c J " u k 1 1 

by Lemma 1431 with D = D l . 
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Next, we turn to the difference of the last two integrals in (14. 191 ). Recalling (14.41) . we write this 
as follows: 



\w ■ nlds, 



V V (p(Fw)P k (w), Jj(tw)= / \w-n\ds- \wn\ds = 
i£J c w eL\{o} JkS ' Js l JdG l 

where n is the outward unit normal to kdVt and dVtk m the first two integrals, while n is outward 
unit normal to dG\. Hence I\{w) > 0, and since G\ C kD l \ D l k , 

0<P k (w)< / \w-n\ds- / \w ■ n\ds < C\w\ 3 , 

where the estimate follows from (14.121) by replacing Vt of Proposition 14. II by D % ; the constant C 
is independent of k. Since the sum X]«,eL\{o} l^^d-^^l) converges absolutely by hypothesis, so 
does the double sum in the previous equation; therefore 



lo {F,n)ds- / 7o (F,n)ds = 0(l). 

dQ.% Jkdn c 

This together with (14.201) and (14.211) shows that R{k) = 0(k 2 ^ 3 ). The normal is irrational a.e. on 

dVl c . Consequently J km 7 <is = J kdn ^ds = k J m ^ds, and (14.161) follows from (14.181) . since 

(1431) holds. □ 

Proofs of the Wulff theorem associated with surface energy minimization llFol over domains 
with given measure typically rely on continuity of the surface energy density with respect to the 
unit normal. Perhaps surprisingly, the extended surface energy density j(F,-): S 1 — > R exhibits a 
dense set of discontinuities as we show next. 

Proposition 4.4. Suppose cp is as in Proposition \3. 3\ and fix F E M^ x2 . Then 

(i) 7 (F, •) : S 1 — > R in (14.41) is Lipschitz continuous on S 1 . 

(ii) 7(-F, •): S 1 — > R defined in (14.151) is continuous at n G S], discontinuous at n G S l R and 
differentiable at most on a subset of S] of measure zero. 

Proof. Arrange the elements of L \ {0} in a sequence: {wj}, j = 1,2,..., such that \wj +1 \ > \vjj\, 
and define gj(n) = (—l/A)Lp{\Fwj\)\wj ■ n\ for n G S 1 . Then clearly gy. S 1 — > R is Lipschitz on 
S 1 and (formally for the moment) 7 (F,n) = YlJLi 9j( n )- Now since \gj\ < Mj = \ip(\Fwj\)\ \wj\ 
on S 1 and the series Y^JLi Mj = J2 w &l\{o} Iv^d-^^DI \ w \ converges in view of Lemma 13.21 then 
Gk{n) = Ylj=i 9j( n ) converge uniformly as k — > oo to j (F, n) on S 1 by the Weierstrass M test. 
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Since n \w ■ n\, n £ S 1 is Lipschitz with constant \w\, the Lipschitz constant of Gk is bounded 
above by 

k 

j2\<p(\ F wj\)\ n < J2 m <o °- 

J'=l u>GL\{0} 

The uniform convergence of the Gk together with the uniform bound on their Lipschitz constants 
guarantee that the limit function j {F, •) is also Lipschitz on S 1 and (i) holds. 



To show (ii), consider the function 



In other words, letting n 




n G Sft {n G M, n/\n\ = n), 
n G S}. 



Then one has 



V2/vi=p/q, {p,q)€Z, gcd(p,g) = l, 
/'('•'i.-'^) = <| o. ^ = 0, ( 4 - 22 ) 

otherwise. 



7 (F, n) = 7o(F, n) + Vn G S 1 . (4.23) 



By (i), it suffices to prove that h is continuous at irrational n and discontinuous at rational n to 
show the continuity part of (ii). In fact, h is very similar to the Thomae function T(x) = 1/q for 
x = p/q, p, q coprime integers (x rational), and zero for x irrational; see e.g. Proposition 4.1 in 
[[Sal . Adapting these results to the h is trivial in view of (|4.22l) . Thus h is continuous at irrational 
n and discontinuous at rational n and so is j(F, ■). Also h is nowhere differentiable by a simple 
adaptation of Proposition 6.1, [[Sal . Since by part (i) 7 D is Lipschitz, it is differentiable a.e. on S 1 
by the Rademacher theorem. Then j(F, ■) fails a.e. to be differentiable by (14.231) . Also it is not 
differentiable at rational n as it is not continuous there. □ 

5 A Continuous Surface Energy Density 

The difficulty caused by the lack of continuity of 7 is in to some extent artificial and can be over- 
come, as we explain next. The previous result suggests that the problem of minimising the integral 
f dn ^(F,n)ds (for fixed F) over a suitable class of domains with fixed may be ill posed, 
in view of llFol (due to lack of continuity of *f(F, •)). However, this is not a physically appropriate 
problem, since by Pick's Theorem, fixing is not equivalent to fixing the total mass, or the number 
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of lattice points of fi. If the minimization were over the class of lattice polygons with fixed lattice 
point number, the appropriate constraint would fix \ + J dQ l/(2\n\)ds instead of Interestingly 
enough, the term l/|n| in the last integral is the one responsible for the lack of continuity of 7. At 
the same time, given a collection of atoms within a convex domain fi, there is some freedom in 
choosing an alternative convex domain f2' containing precisely the same atoms. By choosing in 
a specific way, we can rewrite the surface energy as the integral of a continuous surface density. 
The shortest distance between parallel lattice lines with Miller normal n is l/|n| (also equal to the 
interpalanar distance in 3D). If SI is a lattice polygon, construct Vt' by moving each side with Miller 
normal fii of dtt outward by 1/(2 |n,|), half the inteplanar distance. Then f2' is also a polygon (not a 
lattice polygon) that contains the same atoms as f2 and has the same normal fan, although it is not a 
dilation of f2. Performing the same operation on kVt and writing the energy in terms of the modified 
domain, one arrives at the following representation of the energy: 

Proposition 5.1. Let ip be as in Proposition \3.3\ with d > 3. Suppose Q C M 2 is (a) a lattice 
polygon, or (b) a smooth domain as in Proposition \4.1\ or (c) the piecewise smooth domain of 
Proposition \4.3\ In case (c) assume further that straight and curved sides ofdQ are not tangent at 
their common points. Then for k E Z + , there exists a convex Q(k) C IR 2 containing the same lattice 
points as kft and whose measure equals the lattice point number ofkfl to order 0{k) as k — >■ 00: 

n(k)nL = mm, \n(k)\ = #(knnL)+o(k) = \kn\+o(k), \dn(k)\ = \kdn\+o(i), 

(5.1) 

such that 

E{ktt,y}= [ W{F)dx+ [ j {F,n)ds + o{k). (5.2) 

Moreover, in case (a) the o(k) terms above are 0(1), while in cases (b), (c), they are 0{k 2 ^). 
Finally, in case (b), (15.21) holds for all k 6 R + and not merely integers. 

Proof. Case (a): Suppose VL is a lattice polygon. Then kVt = {x E M 2 : x-n^ < kdi, i = 1, . . . , N}, 
is the intersection of N half-planes of the form x ■ n« < kdi, where n.j e M is the Miller normal of 
the ith side and dj are integers independent of k. Let Cl(k) = {x E M 2 : x ■ fii < kdi + 1/2, i = 
1, . . . , iV}. Thus to construct Q(k), each straight line containing a side of with Miller normal fii 
is translated outward (in the direction fii) by a ^-independent distance l/(2|nj|). The intersection 
of the half-planes of the translated lines is fl(k) This adds to a layer whose thickness equals 
1/(2 1 nil) on the ith side, hence 

|ft(Jfe)| = \m\ + k I 4-ds + 0(1) = #(kVL n L) + 0(1). (5.3) 

J Ml 2 \ n \ 
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The 0(1) term is a correction due to intersection, in the neighborhood of corners, of layers corre- 
sponding to adjacent sides, since directions and thicknesses of layers are /c-independent. The second 
equality above follows from (12.101) of Lemma l2.1[ The 0(1) terms in (15.31) are actually constant (de- 
pend only on Cl and not on k) as is easily shown. This establishes the middle assertion in (15.11) . Since 
the distance of adjacent lattice lines with normal ri, is l/|nj|, the added layers (whose thickness is 
half that distance) contain no new lattice points; thus the first assertion of (15.11) holds true, while the 
last is trivial. Now (15.21) follows immediately from Proposition 13 .31 (15.31) and the definitions (13.41) 
and (l474"l) . 

Case (b): Suppose Q is a smooth domain as in Proposition 14. II Then choose £l(k) = f2, to that 
(157T1) follows from IfHull and note that (1531) is the same as (|43T) with k = r e R+. 

Case (c): Let f2 comply with Proposition l4.3l . For each k let Q(k) be the set obtained by moving 
only the fiat sides kSi C dVtf, i £ Jf of kVl outwards by l/(2|n»|) (and discarding portions of the 
added layers that lie outside the curves Tj near the endpoints where Si join curved sides of <9f2). 
Thus 

\tt(k)\ = \ktt\ +k I -^-ds + 0(1) = #(A;fi n L) + 0(k 2/3 ) (5.4) 

JdQ f 2 N 

The second equality follows from (14.201) and (14.211) . Hence (15.11) holds (the last assertion is easy). 
One again, ([572]) follows easily from (14.171) and (15741 □ 

The second relation in (15.11) means that the measure of Q(k) is chosen to equal the number 
of lattice points of up to the order of the surface energy. This reminds us of the condition in 
Theorem 3 of HBBLII that — in present notation — supposes the existence of a sequence r k — > oo 
such that #(r fe f2 n L) — |r fe f2|. Provided that such a sequence of dilation factors r k exists, then 
it is possible to modify the results of the previous section to show that for the types of domains 
considered here, 

E{r k n,y} = r 2 k [ W(F)dx + r k [ 7o (F, n)ds + 0(1), 
Jn Jan 

thus the non-explicit surface energy density of Theorem 3 in HBBL1 is now determined to be 7 D 
defined in ( 14.41) . Apparently, it does not seem to be known for which domains f2 such a sequence 
of dilation factors r k exists. We are thus led to the construction of the domains Q(k) of Proposition 
15.11 which are not dilations of £1 of the form rfl, since they involve different translations of different 
facets. 

Proposition 15.11 indicates that the appropriate surface energy shape optimization problem of 
Wulff type involves minimizing J dQl j (F, n)ds over a suitable class of domains Q! with fixed. 
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This entails fixing the mass of the crystal. The integrand 7 D is now Lipschitz continuous in the unit 
normal as guarranteed by Proposition l4.4[ 

We must remark, however, that while (15.21) has the aforementioned advantages as regards sur- 
face energy minimization, it is not appropriate as an asymptotic series in k, since the domains of 
integration depend on the latter variable. The appropriate asymptotic series remains (14.161) . 

It is interesting that in cases where the normal is rational on a subset of dQ of positive measure, 
the dilation factors are restricted to be integers. In one dimension it is has been shown HMol that the 
coefficients in the asymptotic expansion of the energy depend on the sequence of dilation factors. 
It should be kept in mind that there is no counterpart in one dimension of an irrational surface, 
which is purely a higher-dimensional occurence. Accordingly, the result of Proposition l4.1l (the case 
of smooth, strictly convex domains with irrational normal almost everywhere on <9f2) is sequence 
independent. 
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